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I. INTRODUCTION 



One of the main subjects of modern physics is to understand the systems consisting 
of many particles interacting with each other. Due to many degrees of freedom, as well 
as due to the interactions between particles, in general, to exactly solve their equations of 
motion is impossible in most cases. However, in some very special cases, that is, for the so- 
called integrable systems, we can find exact solutions for their eigenenergies, wavefunctions 
and other physical properties, because of infinite number of symmetries of the systems. The 
classical examples are dated back to the ansatz, proposed by Bethe in 1931 |I| for the ground 
state wavefunction of the spin 1/2 Heisenberg chain, as well as the Onsager's solution for 
the 2D Ising model in zero magnetic field |@]. 

In 1963, Lieb, and later Lieb and Liniger studied the one dimensional delta potential 
interacting Bose gas wavefunctions by constructing ansatz 0, without knowing Bethe's 
idea for the spin chain. After the work of McGuire on M = 1 case and that of Flicker 
and Lieb on the M = 2 case, Yang, Gaudin independently found the solutions for the 
one dimensional delta function potential spin 1/2 fermion system for M > 2 Together 
with Baxter's solution for the statistical mechanical system (the eight- vertex model) [§, 
factorization of scattering matrix, often called Yang-Baxter relation, has been the key to 
many exactly solvable models, both in condensed matter theory and in high energy physics. 
We'd like to mention, in particular, its great applications to the condensed matter systems, 
such as the well-known Lieb-Wu solution for the one dimensional Hubbard model Pl, the 



Bethe-ansatz for the one impurity Kondo model and the Anderson model 



Having said those, we'd like to mention, as one exciting progress in quantum field theory, 
the development of the conformal field in two dimensions, where the infinite number of 
conformal symmetries fully determine the dynamics of the systems, since it was realized that 
the scaling invariant systems also respect the conformal symmetries ]1^ . The conformal field 
theory has become a powerful tool for the critical phenomenon theory, where the fluctuations 
on all length scales are important and the systems are scale invariant. Generally speaking, 
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such infinite number of symmetries determining systems' dynamics may also be considered 
as one example of solvability when many degrees of freedom are involved. In general, exact 
solvability of many particle systems ( in field theory case, we have continuum of particles) 
have been attracting attentions from both condensed matter physicists, and mathematical 
physicists, as well as from field theorists. 

In this talk, we'd like to view some recent work in the approach to the ID strongly corre- 
lated electron systems by solving them exactly. There has been considerable general interest 
in study of low dimensional strongly correlated electrons, largely due to the discovery of 
high temperature superconductivity ||l3l. The two dimensional t-J model and the Hubbard 
model, beheved to model the high temperature superconductors, have been studied inten- 



sively through various methods [|I4|. The fixed point behaviors of the systems from one 
dimension to two dimensions have been one of the major interests of the theorists working 
in the field. Whether the normal state of the new 2D-like superconducting materials is a 
non-fermi liquid, is one of main focuses for its theoretical study. 

The ordinary one dimensional t-J model with supersymmetry |T5[, and the one dimen- 
sional Hubbard model are exactly solvable [^], with the technique of the Bethe-ansatz. 
Different from the Bethe-integrable systems like these, consisting of one type of integrable 
family, Calogero and Sutherland have developed another quite interesting integrable inter- 
acting many particle family with long range interactions |]30| , pl| . The interests in the 
Calogero-Sutherland type systems are renewed after Haldane and Shastry's independent 
works on the integrable one dimensional spin systems with long range exchange interaction 
in 1988 [l^JT7[] . The close relation of the Calogero-Sutherland type system with quantum 
chaos, random matrices, has been greatly acknowldeged. In our study, however, we focus 
on the integrable correlated electron systems in one dimension related to the Calogero- 
Sutherland- type quantum system, that is, the one dimensional t-J models of hopping 
and exchange. 
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II. ID T-J MODEL ON UNIFORM LATTICE 



Let us first consider the supersymmetric t-J model on a one dimensional lattice of 
equal spaced sites. The Hamiltonian for the system is given by 

N 



H = {1/2)Pg 



J2 Y.^iji^tcja)+ Jij[Pij-i^-ni){l-nj)] 



Pr 



where we take the hopping matrix and the spin exchange interaction to be tij/2 = Jij = 
l/cP{i — j), and d{n) = (L/n) sin(r27r/L) is the chord distance, with L the size of the lattice. 
The operator is the fermionic operator to create an electron with spin component a at site 
i, is the corresponding fermionic annihilation operator. Their anti-commutation relations 
are given by {ci^^,cJ^J+ = {q^,, CjvJ+ = 0, {c^,c]^J+ = 0. We assume that the 

spin component a takes values from 1 to A^. The Gutzwiller projection operator Pq projects 
out all the double or multiple occupancies, Pq = ]\i=i Pci^), and Pci^) = So^m + (^i.n,, with 
ni = Y.a=i4aCia- The operator Pij = E^=i E^=i cLcia'cj^/Cj^ exchanges the spins of the 
electrons at sites i and j, if both sites are occupied, and nj are the electron number 
operators at sites i and j. 

The SU{2) model was first introduced by Kuramoto and Yokoyama, who provided the 
ground state of the system, the Gutzwiller projected Fermi sea away from half- filling ||21|| . 
More generalized Jastrow wavef unctions, including the ground state as a member, were also 
introduced The asymptotic energy spectrum was derived by Kawakami based on 

the assumption of asymptotic scattering matrix factorization ^2|. The system is identified 



as a free Luttinger hquid [I21| , p2|,|26| , p3| . In particular, for the SU{2) case, the asymptotic 
Bethe-ansatz spectrum was explicitly shown to be exact, and the correct thermodynamics 
was given when the spinon rotation was properly taken into account [^], in terms of a set of 
generalized Laughlin- Jastrow wavefunctions. In general, exact solvability implies existence 
of infinite number of constants of motion. At half-filling, this system reduces to the Haldane- 



Shastry spin. Quite recently, the constants of motion of the spin chain have been studied [28 



using the operators Mij . In our most recent work, we have pointed out that their results 
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can apply to the case where the holes are involved, with the superalgebra representation of 
the fermionic fields in the t- J model |^ , by mapping the system to a mixture gas of fermions 
and bosons on the one dimensional lattice. 



A. The constants of motion 

For the t-J model, we may introduce two new fields, the / and b fields. For the new 
fields, we have {fia,fja'}+ = 0,{/t,/jv}+ = 5ij5^^> ,[bi,bj] = 0, [foi,^^] = 5ij. The b field 
always commute with the / field. The size of the Hilbert space at each site is oo in this case. 
After projecting out the zero occupancy, double and multiple occupancies, the new Hilbert 
space is equivalent to the one defined by the c fields with no double or multiple occupancies. 
In particular, we may represent the fermionic electron operators cl^ and Qo- as: 

PG{i)Ci„PG{i) = '^l,nj,+n}^l/i<T'^l,n'j+n}) (2) 

where + = + J2a=i fLfio-- Hence, a state vector can be written as 
10 >= J2 <Pixiai,X2a2,- ■ ■ ,XN,aNjyi,y2,- ■ ■ ,yQ) X 

CTi,(72,---,CTjVg {x},{y} 

Xfliai fL<T2 ' ' ' f^N, ajve ^yi ■ ■ ■ I > ' (^) 

where is the number of / fermions on the lattice, Q is the number of b bosons, and we 
require that Xi 7^ xj Vk 7^ Vi, and that the / fermion positions {x} and the b boson positions 
{y} span the whole chain. Obviously, Ne is also the number of the electrons, and Q is also the 
number of holes on the lattice. The amplitude is anti-symmetric when exchanging (xjCij) 
and (xjaj), and symmetric in the boson coordinates {y} = {yi, y2, - ■ ■ , yq). The Hamiltonian 
of the supersymmetric t-J model, using the mapping Eq. (|^) in a straightforward way, can 
be written in terms of the fermionic / field and the bosonic b field. 

With the above mapping, we may write the eigen-energy equation of the supersymmetric 
t-J model in first quantized form. Following Ref. ( we define the "exchange opera- 

tor" Mij as MijF{qi, q2, ■ ■ ■ ,qi, ■ ■ ■ ■ ■ ■ , Ql) = F{qi, q2, ■ ■ ■ ■ ■ ■ ■ ■ ■ , Ql), where the 



function F is an arbitrary function of some position variables (gi, q2, ■ ■ ■ , Ql), i-e., the oper- 
ator Mij exchanges the positions qi,qj of the particles i and j. In terms of such exchange 
operators, the eigen-energy equation of the t-J model takes the form as follows 



(1/2) 



E d-\q,-qj)M,, 

l<i¥^j<L 



0({g};{a}) = E0({g};M), (4) 



where {q} = (gi, ^2, ■ " " , 9l) = (xi, Xs, ■ ■ ■ , x^v,, I/2, ■ " " , and {ex}) = 

4>{(ii'^i,(i2<^2,- ■ ■ ,qN,<^NjqN,+iqNe+2- ■ -ql) = 4>{xi(Ti,x2(j2, - ■ ■ ,XN,(jNjyi,y2, - ■ ■ ,yQ) is the 

amplitude of the state vector of Eq. (^). {a} = {cri,o-2, ■ ■ ■ ,o'nJ are the spin vari- 
ables of the / fermions. The operation Mij is defined in the conventional way: 
Mij(l)i{q};{cr}) = {o-}), with {q} = {qi, q2, ■ ■ ■ , Qu ■ ■ ■ , Qj, ■ ■ ■ , Ql) and {q'} = 

ilij Q2, ■ ' ' yQjr ' ' yQir ' ' : Therefore, the original t-J model is written as an eigenvalue 
problem for a mixture of the / fermions and the spinless h bosons in terms of the position 
exchange operator Mij. 

The constants of motion have been studied recently for the spin chain of exchange 
interaction. In the work, a system of identical bosons on the chain is considered, and 
the wavefunctions is totally symmetric when interchanging two particles. In terms of the 
operator Mjj, the following commutation results have been found 

[/n,/J=0, (5) 

where /„ = Ef=i ^rf , with tt^ = Y.j(^i){zj / Zij)Mij, Zi = e^"""^'^^ , Zij = Zi - Zj, and n,m = 
0, 1, 2, ■ ■ ■ , 00. It was found that all these quantities commute with the Hamiltonian H = 
J2i<ij^j<L l^i ~ Zjl^'^^ij as long as the particles occupy the whole chain. For a system of 
identical bosons on the chain, the wavefunction is totally symmetric when we simultaneously 
interchange spins and positions of two particles. The effect of the exchange operator Mij 
is just equivalent to the effect of the spin exchange operator alone. With this method, the 
constants of motion for the SU{N) Haldane-Shastry spin chains were constructed in Ref. 



In the recent work El], we have stressed that, in the language of the exchange operators 



Mij, the commutation results proved by Fowler and Minahan hold as operator identities. 
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The central issue is that the forms of wavefunctions of many particle systems, as well as 
the statistics of the particles or the types of the particles, do not matter in order for the 
commutators to hold, as long as the particles occupy the whole chain. We may then apply 
the "exchange operator formalism" to the wavefunctions of mixtures of fermions and bosons. 
Therefore, from the eigen-equation Eq. (H), we thus conclude that in the first quantization, 



all the invariants of the t-J model are the same J^'s as in Ref. 28 



With the permutation properties of the amplitude for the mixture of bosons and 
fermions, it is straightforward to write all the invariants of motion of the t-J model in the 
second quantization form using the /„'s. For instance, the exchange operation between the 
/ fermion positions is equivalent to the spin exchange operation ( minus sign involved ), 
the exchange operation between h boson positions is equivalent to the hole-hole interaction 
term, and exchange operation between / fermion and h boson positions is equivalent to the 
electron hopping. Such procedure to reduce an J„ to a second quantized form is quite simple, 
and we do not write all the details. Thus we provide a systematical way to construct all the 
invariants of motion for the super symmetric t-J model, either in first quantized or in 
second quantized forms. 



B. The Jastrow product wavefunctions 

Having said the integrability of the t-J model, we'd like to mention that the exact eigen- 
wavefunctions can be constructed explicitly for the system. In the SU{2) case, states in 
the Hilbert space can be represented by spin and hole excitations from the fully-polarized 
up-spin state \P >■ If we let Q denote the number of holes and M denote the number of 
down-spins, then Sz is given by Sz = {N — Q)/2 — M. The wavefunctions are given by 

m=T.i^ix,y)Y[S;^l[hl\P), (6) 

x,y a i 

where the amplitude 'ip{x, y) is symmetric in x = (xi, X2, . . . , xm)^ the positions of the down- 
spins, and antisymmetric in y = {yi,y2, ■ ■ ■ ,yQ), the positions of the holes. = c^.^|C^^-|- 
is the spin- lowering operator at site Xa and /ij. = Cy-^ creates a hole at site 
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We can construct a general class of states corresponding to states of uniform motion 
and spin polarization, by generalizing Kuramoto and Yokoyama's Jastrow ground-state as 
follows 



'^Gi.x.y; Js, Jh) 



exp 



*o(a;,y) = n 



*o(a;,y) 

a Vi] 



(7) 



a</3 i<j 

Here, Jg and Jh govern the (uniform) momenta of down-spins and holes respectively. Jg and 
Jh take on either integral or half-integral values as appropriate to insure that tpc has the 
correct periodicities under Xa Xa + N and Ui ^ yi + N . 



The wavefunction tpc is found to be an exact eigenstate of H with energy 



/V2 2 

—E = -M{M^ - 1) - 2MJs{N - Js 



3 
+Q 



i(iV2 - 1) + ^(Q2 - 1) + i(M + Q)(2M - Q) 



-2J;,(iV - Jh) + 2{Js - Jh 



+ ^{N'-l){N-Q). 



(8) 



which is valid under the conditions | - N/2\ < N/2 - (M - 1 + Q/2), \Jh - N/2\ < 
N/2~{M + Q — 1)/2 and \ Jh — Js\ < (M + l)/2. For a given 5*;^, the minimum energy is given 
when Js and are as close to N/2 as possible. When Q = this reduces to the result for the 
Heisenberg chain [0,|T^. These energy levels were also found independently in Ref. p3[. The 
equivalence of the two representations is clear when identifying = Jh — N + {M + Q + l)/2 
and Ji = Jh — Js — {M — l)/2. Besides these explicitly constructed excited wavef unctions, 
we'd like to emphasize that the full spectrum and the thermodynamics can also be written 
out in terms of Jastrow wavefunctions of more generalized form |^ . 



III. THE ID T-J MODEL ON NONUNIFORM LATTICE 



A. The integrability 



Besides the above integrable t-J model on equally-spaced sites, let us consider another 
supersymmetric t-J model of hopping and exchange on a chain with sites not equally 
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spaced. The positions of the sites xi^X2, ■ ■ ■ are determined by the equation 

Xi= 'i/{x^-Xjf. (9) 

An integrable spin chain was introduced on a non-uniform lattice on which the sites are 
positioned by the equation Later, Frahm reahzed that the solutions of the equations 



are the zeros of the Hermit polynomials Hl{x) |2^. The Hermite polynomial is defined as 

H,ix) = (-l)V^^^- (10) 

It was well-known that the Hermit polynomial has L roots, all of which are real numbers, 
and the lattice is thus well defined. 



Doping the spin chain, we have introduced the following supersymmetric t-J model p3 

N 



H = (1/2)Pg 



Pg, (11) 



where the hopping matrix and the anti-ferromagnetic exchange interaction are given by 
tij/2 = Jij = l/{xi — XjY, and each site is occupied at most by one electron. 

In the half filling case iVg = L, this system reduces to the spin chain that has been studied 
before, for which a similar exchange operator formalism has been developed. Let us just write 
down the results obtained in Ref. [^]: [/„,/,„] = 0, [IniH] = 0, where /„ = Y.i=i ^« = 
a|aj, and a| = nj+iqi, ai = 7ii-iqi, with 7ii = T.j(^i) i{qi-qj)~^Mij, H = Y.ijLj{qi-qj)~'^Mij, 
and n,m = 0,l,2,---,cx). Here, all the particles are put on the chain where the sites are 
positioned as determined by the Eq. @. We may relate the operation of exchanging particle 
positions to the operation to exchange particle spins, by assuming that we have identical 
bosons again, for which the wavefunctions are totally symmetric when we exchange the spins 
and positions of two particles simultaneously. With this assumption, from J^'s, we thus can 
derive all the invariants of motion written in terms of the spin exchange operators alone. 

Again, all commutation results written in terms of the exchange operators Mij obtained 
by Polychronakos hold as operator identities, as long as the particles occupy the whole chain 
of the sites positioned in the special way. The forms of the wavefunctions do not matter at all. 
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Thus the commutation results may apply to wavefunctions of particles of arbitrary statistics 
or wavefunctions of mixtures of particles of different statistics on the chain. Mapping our 
new supersymmetric t-J model in terms of the b and / fields, we may also write the eigen- 
energy equation in first quantized form. In terms of the exchange operators between the 
positions of the bosons and fermions, the Hamiltonian takes the form 

if =(-1/2) ^ (g,-g,)-2M,,. (12) 

Applying the formalism to this t-J model, in a similar way we may obtain all the invariants, 
either in first quantized or in the second quantized form, which commute among themselves 
and with the Hamiltonian. Thus this supersymmetric t-J model is also completely integrable. 



B. The wavefunctions 

The Hamiltonian commutes with the permutation operator Tij = P"-Mij exchanging the 
/ fermion spin and position simultaneously. Let us work in the Hilbert space where the 
number of fermions of each fiavor is fixed, i.e. N^j, o = 1,2, . . . , N, is fixed. Gruber and 



Wang have considered the following wavefunction in Jastrow product form 



where {x} and {y} span the whole lattice. 



It was shown that the wavefunction is an eigenstate with eigenvalue [25| 



N 

E = -L{L - l)/4 + (1/2) Y^iN^ - 1)N^. (14) 

(T=l 

Although it is expected that this wavefunction is the lowest energy state in the subspace of 
fixed Ni, N2, ■ ■ ■ , N^, we were not able to prove it. However, in the case of SU (2), the small 
lattice diagonalization up to 8 sites confirms this conclusion. 

For fixed number of the electron number, the minimum of the energy is obtained 
when lA^^o- — A^o-'l is as small as possible for each pair a 7^ a'. In the SU{2) case, the above 
result becomes 
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E = (-l)L(L - l)/4 + {1/2)N^{N^ - 1) + il/2)N^{N^ - 1), (15) 

where N-^ and iV| are the numbers of the up-spin electrons and the down-spin electrons 
respectively. For fixed number of electrons on the chain, i.e. for fixed A^e, the minimum of 
the energy given in Eq. (|l5|) is obtained when Sz = for even A^e, or when 5*2 = ±1/2 for 
odd Ng. Therefore, the ground state is a spin singlet ( respectively spin 1/2 ) state for even 
( respectively odd) number of electrons on the chain. In particular, for an even number of 
electrons on the chain, the ground state energy is 

Eg = (-1/4)L(L - 1) + iVf/4 - NJ2, (16) 

while for an odd number of electrons it is 

Ea = {-l/mL-l) + {^^r. (17) 

The charge susceptibility of the ground state Xc is given by Xc^ = d'^Ec/dN^ = 1/2, 
independent of the electron concentration. The charge susceptibility is also finite at half- 
filling Ne = L, in spite of the existence of a metal-insulator phase transition at half-filling 
for this system. 

To study the spectrum of the system away from half-filling, we follow the idea introduced 



in Ref. ^7^, defining the operators: 



= ^ II {qj - Qk) ^Mjk = tt], 
Oj = 71 j + iqj, 

ttj = 71 J - iqj, (18) 
which satisfy the following commutation relations: 

[7Tj,7lk] = 

[qj,H] = i7ij 

[n„H] = -2z ^(g^.-g,)-3. (19) 
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Then using the property of the roots of the Hermite polynomial we have 

K-,^] = -%-, (20) 

and 

[a], H] = —a] 

[aj,H]=a,. (21) 

Therefore the operators AKu) = a] Sl''\i = 1, 2, ■ ■ ■ , iVe, where u = 0,±,z for the SU{2) 
case with S^^^ = 1, will act as raising operators, while their hermitian conjugate v4j(z/) will 
act as lowering operators. It thus follows that the wavefunction 

'^W,M = En(4(^.))"^'0 (22) 

P i=l 

with {n} = {rii, n2, ■ ■ ■ , nNj,ni > 0, {u} = (z/i, z/2, ■ ■ ■ , i^nJ, is either an eigenstate with 
energy 

E{n} = E + Y,ni (23) 

i=l 



or zero. Moreover, for the b boson degrees of freedom, we have following property 



a,(/) = 0,ze (A^e + l,iVe + 2,---,L), (24) 



We have also shown the following results are true: 



L 

aa(j) = 0, 

a=Ne+l 

Y.ai(j) = 0. (25) 
1=1 



Furthermore, we realize that 



iVe 

[Y^A{z)]^ = 0. (26) 

i=l 

In the particular case where N-^ = Ni, the wavefunction function is a spin singlet and we 



may globally rotate Eq. (^) in the spin space, giving us p5 
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i=l 



0. (27) 



These results have demonstrated the impossibihties of constructing non-vanishing eigen- 
states with these lowering operators and the wavefunctions (f) in the subspace where the 
number of electrons of each flavor is fixed. The results are in support of our conclusion 
that in general the Jastrow wavefunction is the lowest energy state in the subspace of fixed 
number of particles of each flavor. 

In the end, we have concluded that the excitation spectrum of the system is of the form 



Pi 



E{s) =E + s- 

s e (0,1,2, (28) 

i. e., the spectrum of this t-J model consists of equal-spaced energy levels. Since the model is 
on a finite chain, Smax is finite. In the special case of SU{2), the small lattice diagonalization 
up to 8 sites suggests that the highest energy level is given by 

EmaxiQ) = L{L - l)/4 - Q{Q - l)/2, (29) 

where Q is the number of holes on the chain. The feature of the t-J model spectrum consisting 
of equal-distant energy levels may also be seen by taking the strong interaction limit of the 
Calogero-Sutherland-Morse quantum system for a mixture of fermions and bosons, 

H = (-1/2) j: d'/dq^ + E + E ^(^ - M.,)/(g. - g,)^ (30) 

i=l i=l i<j 

where there are Ne fermions with spins and Q spinless bosons, Mij permutes the positions 
of the particles i and j only. 

Finally, we'd like to point out that the Hilbert space of the t-J model has very interesting 
degeneracy pattern, as indicated by the finite lattice study. It still remains to find a system- 
atic rule to characterize the degeneracy of each energy level. So does the thermodynamics 
remain to be solved. 
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IV. CONCLUDING REMARKS 



In this talk, we have reviewed some recent progresses in study of the ID strongly corre- 
lated electron systems. They are exactly solvable systems, close related to the continuum 
Calogero-Sutherland quantum nonrelativistic system with long range pairwise interaction. 
There have still been constant interests in such systems recently. In particular, we would 
like to mention that the Calogero-Sutherland system can be interpreted, with Haldane's 
new definition of statistics based on counting state number, as a system of particles obeying 
fractional statistics, where the dynamic coupling is closely related to the fractional statisti- 



cal parameter |l36|-pq]. We also wish to mention that the complete analytical confirmation 
for the energy spectrum and thermodynamics of the 1/r Hubbard model with finite U, 
proposed by Gebhard and Ruckenstein have been still lacking |^. In the strong coupling 
limit U = oo, the Gutzwiller-Jastrow wavefunctions were demonstrated to be eigenfunctions 
of the Hamintonian, the spectrum conjectured by Gebhard and Ruckensten was explicitly 



proved to be exact ||3^ , |35[| . It remains to find the wavefunctions and the constants of motion 
for the finite U case for this Hubbard model with only relativistic right movers. Further 
exciting results are expected for these interesting integrable systems from the communities 
of condensed matter theory, high energy physics, as well as of mathematical physics. 
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